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Abstract 

Recently  developed  formulations  of  the  inverse  problems 
in  acoustics  and  electromagnetics  are  described.  There  are  two 
types  of  formulations,  one  in  the  geometrical  optics  limit  and  the 
other,  an  exact  formulation  for  the  inverse  source  problem.  Both 
basic  formulations  are  extended  to  include  the  realistic  problem 
of  a "limited  aperture"  of  observations.  It  is  also  shown  that 
the  inverse  source  formulation  can  be  applied  to  the  problem  of 
reconstruction  of  media  inhomogeneities  from  remotely  sensed  field 
data.  The  basic  physical  optics  result  is  that  the  characteristic 
function  of  the  scattering  obstacle  and  the  phase  and  range  normal- 
ized scattering  amplitude  are  a Fourier  transform  pair.  All  other 
formulations  lead  to  Fredholm  integral  equations  of  the  first  kind. 


INTRODUCTION 


The  objective  of  inverse  scattering  methods  is  to  obtain 
information  about  a source  distribution,  a scattering  obstacle 
or  inhomogeneity  of  a medium  from  remotely  sensed  field  data. 
Recently,  we  have  formulated  such  inverse  problems  in  acoustics 
and  electromagnetics  in  terms  of  Fredholm  integral  equations  of 
the  First  Kind.  Furthermore,  some  of  these  formulations  take 
account  of  the  realistic  difficulty  that  the  f Ids  are  usually 
observed  only  over  some  " limited  aperture",  beta  spatially  and 
temporally,  rather  than  in  all  directions  and  for  all  time. 


In  this  paper  we  shall  describe  these  formulations  and 
present  the  available  evidence  of  their  validity.  This  evidence 
is  admittedly  sparse  at  this  time,  but  is  also  extremely  promising. 
Our  purpose  is  to  collect  these  recent  results  and  to  communicate 
them  because  we  believe  there  is  a wealth  of  research  directions 
to  be  explored  well  beyond  our  own  capacities  and  interests. 


The  list  of  applications,  non-military  as  well  ac 
military,  is  quite  extensive.  Among  these  are  the  following. 


(1) 

(2) 


Construction  of  images  of  tumors. 


Analysis  of  subsurface  strata  for  resource 
identification  and  recovery. 


(3) 


Location  and  identification  of  discharges  in  storms 
for  the  analysis  of  the  storms  themselves  and 
prediction  of  tornadoes  from  charasteristic  source 
patterns. 


(4)  Location  of  buried  bodies — an  important  problem 
in  law  enforcement. 


m 
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(5)  Construction  of  images  of  airplanes,  missiles, 
surface  vessels  and  submarines. 


(6) 

(7) 


Landmine  detection. 


Analysis  of  soundspeed  variation  (and  consequent 
sonar  degradation)  in  the  water  and  sea  bed 
sedimentary  layer. 


We  shall  see  below  that  the  method  of  deriving  the 
integral  equations  is  not  peculiar  to  the  wave  equation  or 
Maxwell's  equations.  Thus,  it  is  likely  that  for  wave  phenomena 
governed  by  other  equations,  completely  analogous  formulations 
of  the  inverse  problem  are  possible. 


There  are  two  distinct  categories  of  problem  formulation 
to  discuss,  a physical  optics  (approximate)  formulation  and  an 
exact  integral  equation  formulation.  The  former  evolved  from  a 
basic  identity  derived  a number  of  years  ago  by  the  second  author 
(Bojarski,  1967,  and  Lewis,  1969)  for  the  time-harmonic  acoustic 
case.  The  latter  evolved  from  more  recent  work  of  his  (Bojarski, 
1973)  for  time  harmonic  wave  propagation. 


The  physical  optics  identity  relates  the  range  and  phase 
normalized  back  scattered  (monostatic)  far-field  cross-section 
to  the  Fourier  transform  of  the  characteristic  function^  of  the 
scattering  body.  This  result  has  now  been  extended  to  the 
bistatic  case.  In  this  form,  much  more  information  is  extracted 
with  only  two  stationary  transmitters  and  an  array  of  receivers. 
These  results  are  described  in  Section  II  and  a nunerically 
synthesized  check  on  the  method  reported  by  the  second  author 
(Bojarski,  April  1973,  February  1974)  are  presented  in  Section  IV. 


In  another  extension  discussed  i.-  Section  III,  we 
dispense  with  the  far-field  approximations  and  obtain  an  integral 
equation  for  the  characteristic  function.  This  formulation  would 


TThe  characteristic  function  is  unity  inside  the  scatterer  and 
zero  outside,  thus  describing  the  scattering  body. 
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be  appropriate  when  the  distance  between  the  observed  remote  field 
and  the  obstacle  was  comparable  to  the  length  scale  of  the  obstacle. 

The  consequences  of  limited  aperature  for  far-field  physical 
optics  are  discussed  in  Section  IV.  Again,  numerical  results  by 
Bojarski  (April  1973)  are  described. 

In  the  exact  formulation,  one  obtains  an  inhomogeneous 
integral  equation  for  a source  distribution  with  the  inhomogeneous 
term  being  an  integral  of  the  data  observed  remotely.  Bo jar ski's 
basic  integral  equation  for  wave  propagation  governed  by  the  Hemoltz 
equation  is  described  in  Section  V.  Again,  numerically  synthesized 
checks  on  the  basic  formulation  are  presented.  Significant  among 
these  is  an  example  in  which  two  point  sources,  one-half  wave  length 
apart  are  resolved. 

Again,  the  limited  aperature  problem  is  of  great  interest 
and  the  formulation  which  takes  account  of  this  feature  is  presented 
in  Section  VI. 

Problems  in  which  one  seeks  an  inhomogeneity  of  a medium 
can  be  reformulated  as  inverse  source  problems  amenable  to  analysis 
by  this  exact  formulation . Inverse  scattering  problems  for  "trans- 
lucent" obstacles  (neither ^acoustically  hard  nor  soft)  can  also  be 
viewed  as  inhomogeneity  problems  and  are  thus  reduced  to  inverse 
source  problems.  The  "source"  here  is  not  to  be  confused  with 
"equivalent  sources"  which  are  physically  unrealizable  mathematical 
idealizations.  These  ideas  are  expanded  upon  and  clarified  in  Sec- 


tion VII. 


Extensions  of  the  inverse  source  problem  formulation  to  the 
time  domain  are  discussed  in  Section  VIII.  These  time  domain  formu- 
lations the  inherent  advantages  noted  above.  Furthermore,  there 

are  two  distinct  formulations  depending  on  whether  or  not  the  source 
distribution  decays  to  zero  sufficiently  rapidly  in  time.  When  this 
is  not  the  case,  it  is  seen  that  the  time  harmonic  case  is  not  neces- 
sarily equivalent  to  a transform  of  the  time  dependent  formulation. 

Ml  of  the  integral  equations  obtained  here  are  Fredholm 
intf?* •u  .quations  cf  the  First  Kind.  The  problems  are  mathema- 
tica_jL^  ill-posed,  suffering  both  from  non-uniqueness  and  non- 
continuous  dependence  on  the  data. 

As  to  the  former  difficulty,  we  believe  that  the  non- 
uniqueness resides  in  a class  of  non-phvs.ically  realizable  solutions 
such  as  source  distributions  "at  infinity"  or  scattered  waves  con- 
taining uigenmodes  of  annular  regions,  which  modes  are  not  "outgoing". 
This  is  a conjecture  to  be  verified  or  contradicted. 

The  ill-posedness  has  as  practical  implication,  the  fact 
that  the  solutions  are  "noise-limited".  That  is,  relatively  small 
amounts  of  noise  can  degrade  the  solution  so  much  that  the  "true" 
solution  is  unrecognizable.  The  quantification  of  this  observation 
is  well  worth  studying. 

The  extension  of  the  formulations  presented  here  to  Max- 
well's equations  is  straightforward.  They  can  be  found  in  the  above 
cited  Bojarski  papers.  Since  the  objective  of  this  report  is  to 
collect  the  basic  formulations,  we  do  not  repeat  those  extensions 
here. 


with  ug  the  scattered  field.  It  is  required  that  uT  satisfy  the 
equation 


(V2  + k2)  uT  = 0 


(2.3) 


and  the  boundary  condition 


uT  = 0 on  SB 


(2.4) 


The  scattered  field  must  satisfy  the  same  equation  (2.3), 
must  be  such  that  (2.4)  is  satisfied  and  be  "outgoing”  at  infinity. 
One  can  then  derive  the  following  integral  equation  for  Ug 
(Baker  and  Copson,  1950). 


=/  k<2’'S>  f?fS‘5 

J 3BT 


3tu  (x*) 

g<*-K’)-s#  " 


ds1 


(2.5) 


Here  SB  denotes  the  boundary  of  B,  S/Sn1  denotes  the 
normal  derivative  and 


x)cx 

g{x)  = ? x = |x|  (2.6) 


is  the  free  space  outgoing  Green's  function. 

The  physical  optics  approximation  consists  of  making  the 
following  assumptions  about  the  fields.  The  incident  plane  wave 
illuminates  a portion  of  SB  called  the  lit  side  and  denoted  by  L 
in  Figure  2.1.  We  may  characterize  L analytically  by  the 


condition  k*fi>  o cn  L,  with  fi  the  outward  unit  normal  to  B.  The  remainder 
of  3B  is  denoted  by  D in  that  Figure.  It  is  then  assumed  that  Ug  satifies 
the  following  approximate  boundary  conditions 


Ug  = - -exp  {i  k*xf 


3uS 

3n 


3uI 

W 


= -ik  • fi  exp  {ik*xf 


on  L , Uo  = ^ = 0,  on  D 

® an 

(2.7) 


We  may  then  substitute  (2.7)  into  (2.5)  to  obtain 


Ug  (x,k)  = ij  exp  | ik*x'  | | (x-x' ) + ik*fi  g (x-x ' ) j>  dS  (2.8) 


wbsre  we  have  used  the  explicit  form  (2.1)  for  Ug. 

We  now  proceed  to  introduce  the  far- field  approximations 
into  this  result.  Firstly,  we  set 


Ug(x,k)  = p(k)  exp  (i  kx{  /4ttx.  (2.9) 

Here,  the  function  p(k)  is  known  as  the  range-and-phase  normalized 
far-field  anplitude. 

Secondly,  we  expand  |x-x'|  for  x»x'  and  use  this  result 

in  g: 

g (x-x' ) = exp  {ik  [x-ik.x'Jj  / ^ ; 


a 


g (x-x')  = ik&*fi  exp  { ik  [x-iibx'jj  /4mx 


(2.10) 
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Here  ft  is  a unit  vector  in  the  direction  of  x. 


The  results  (2.9)  and  (2.10)  are  substituted  into  (2.0) 
and  the  common  factor  exp  {ikx}/4Trx  is  removed  to  yield 


p(k,ft)*  - JiH-  £k  - kft]  exp  {ix' • |k  - kft]}  dS' 


(2.11) 


Let  us  now  illuminate  the  scatterer  with  the  incident 


wave, 


Uj  (x,  - k ) = exp  {-i  k • x}. 


(2.12) 


and  observe  at  the  point  -x.  The  range-and-phase  normalized 
scattering  amplitude  for  this  case  is  then  p (-k,  -ft) . In  this 
case,  ic  is  D that  is  illuminated  and  p<-*k,-ft)  has  a representation 
similar  to  (2.1)  except  that  the  k and  ft  have  changed  sign: 


p(-k,  -ft)  = j ifi  • [k  - kft]  exp  {-ix'  • [k  - kft]  } dS* 


(2.13) 


We  take  complex  conjugate  (denoted  by  *)  in  this  equation  and 
add  to  (2.11)  to  obtain 


p(k,ft)  + p* (-k,  -ft)  = - J ifi  • |k  - kft]  exp  fix-  £k  - kft]}  dS' 


(2.14) 


L & D 

The  integrand  here  can  be  continued  to  the  interior  of 
B as  a regular  function.  Thus  we  may  apply  the  divergence 


theorem  to  rewrite  the  right  side  as  a volume  integral: 


p(k ,St)  + p*(-k,  St)  = jk-kfc]2  j exp  {ix*.  j}c— k5tj  }dV! 


{2.15) 


Let  us  denote  by  Y(x)  the  characertistic  function  of 
the  scatterer  B and  simply  rewrite  (2.15)  as 


p(k,x 


= /*y  (x1 ) exp  {ix’ * [k-k&l  } dV* 
(k-kx)2  I - L-  J 


(2.16) 


where  now  the  domain  of  integration  is  all  of  space. 


We  observe  from  (2.11)  and  (2.13)  that  both  o(k,  St) 
and  p*(-k,  -St)  are  functions  of  k -k&.  Thus,  let  us  define 


k = k-k& 


(2.1.7) 


p(k,ft)  + p*  (-k,-&) 


(2.18) 


Then,  (2.16)  can  be  rewritten  as 


Iy~ 


?(k)  = lY(x')  exp  {ix'*K}  dV* 


(2.19) 


We  see  here  that  the  observed  fields  have  produced  the 
Fourier  transform  of  the  characteristic  function  of  the  scatterer 


If  Y(£)  is  known  for  all  £,  then  we  can  invert  (2.19)  to  yield 


Y(x)  = -^yj  J 9('k)  exp  {ix**}  dx3 


In  the  earlier  back  scattering  derivation. 


kk  - -k;  k - 2k 


Thus,  for  this  case  one  illuminates  and  receives  at  all  wave 
numbers  and  in  all  directions  in  order  to  cover  all  k . For 


the  general  result  in  which  jc  is  given  by  (2.1?)  rather  than 
(2.21),  the  obstacle  is  illuminated  from  only  two  directions, 
-k  and  at  one  wave  number  k.  Then,  observations  in  all 
directions  fill  out  a sphere  in  £ space  with  center  at  k and 
radius  k,  as  shown  in  Figure  2.2. 


FIGURE  2.  2 


It  is  well  known  in  radar  analysis  that  the  physical 
optics  approximation  is  poor  in  the  bistatic  case  if  the 
directions  of  transmission  and  observation  are  too  widely 
separated;  i.e.,  if  the  angle  between  k and  -k  becomes  too 


~ i Tr  1 ••a 


large.  Nonetheless#  for  a cone  of  small  angle,  such  as  is 
shown  in  Figure  2.3,  the  bistatic  result  (2.17,  19)  still 
remains  valid.  In  conjunction  with  the  limited  aperture 
discussion  in  Sect^onlV,  this  observation  should  prove  quite 
useful. 


FIGURE  2.  3 


11 


In  this  section  we  shall  develop  an  integral  equation 
for  the  characteristic  function  of  a scattering  obstacle  with 
the  aid  of  the  physical  optics  approximation  but  under  the 
constraint  that  the  fields  are  observed  too  near  the  scatterer 
for  the  far-field  approximations  to  be  valid. 


FIGURE  3.  1 

We  suppose  that  the  scattered  field  and  ics  normal 
derivative  are  observed  on  the  boundary  3£>  of  the  domain  f)  of 
Figure  3.1.  As  in  the  previous  section,  we  assume  that  is 
the  response  to  the  incident  field  Uj  of  (2.1)  subject  to  the 
conditions  (2.2,  3,  4).  The  result  (2.5)  is  still  valid  as 
well*  Now,  however,  we  also  write  down  the  solution  representation 


“Ug (x,k)  = /{ Ug(x'k)  (x-x' ) - g* (x-x' ) (xjk)  dS  + 0(x,k) 

J ah*  to' 


(3.1) 


8B 


x e f> 


Here  (*)  denotes  complex  conjugate  and 
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0{x,k)  = - j{ 


{vu  (x*  ,k' ) j&i.  (x-x' ) - g*  (x-x‘ ) 
® “ 9n* 


aus 

3n‘ 


<x‘  ,k) } dS' 


3p 


(3.2! 


The  integral  0 arises  in  this  representation  because  g*  does 
not  satisfy  the  radiation  condition  that  g does.  We  note, 
however,  that  0 contains  in  it  the  data  observed  on  3p  and  is 
thus  a known  quantity. 

By  subtracting  (3.1)  from  (2.5)  we  obtain 


f 


(t^(x',k)  (x-x')  - gjOc-x’) 


aug 


3n' 


(x'  ,k) } dS'  = 0 (x,k) , x e p 


(3.3 


3B 

Here 


g-j-fx)  = g(x)  - g*(x)  = —S"*  * 


(3.4 


We  now  make  the  physical  optics  approximations  (2.7) 
and  obtain 


j ft  • V |gJ  (x-x' ) exp  (ik*x* } j dS'  = G (x,k) , x 


£ P 


(3.5 


As  in  Section  II  we  now  carry  out  the  same  analysis  for 
th^  incident  wave 


UI  (x,“k)  = uT*  (x,k)  = exp  {-i  k-x} 


(3.6 


-■*>< 


- Jk'V  |gj (x-x1 ) exp  {-i  k*x}|  dS'  = 0(x#-k),  xe  D 


(3.7) 


the  domain  of  integration  being  D here  because  this  is  the 
surface  illuminated  by  the  incident  field  (3.6). 

We  take  complex  conjugate  in  (3.7) , using  (3.4),  and 
subtract  this  result  from  (3.5)  to  obtain 


w 


I 


^fvV  jgjtx-x')  exp  {ik*x'}j  dS'  = Oj(x#  k ) , x ep 


(3.8) 


9B 


Here 


0j(x,k)  = 0 *(x,-k)  - 0(x,k) 


(3.9) 


ft 


m 


As  in  Section  II,  we  can  now  apply  the  divergence 
theorem  to  rewrite  the  left  side  of  (3.8)  as  a volume  integral. 
Indeed,  here  and  in  the  previous  section,  it  is  the  physical 
optics  approximation  which  makes  this  crucial  step  straight- 
forward. The  result  is 


j |gj(x“x' ) exp  (ik  x'}|  dV'  = O^Xfk),  xej) 


fi 


(3.10) 


We  may  rewrite  this  equation  in  the  form 
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(3.11) 


jy2  |gj (x-x' ) exp  Uk*x' }j  y(x')  dv'  = 0j(x,k),  x e V 


This  is  a Fredholm  integral  equation  of  the  First  Kind, 
the  first  of  many  to  be  encountered  here.  It  is  well  known 
that  these  equations  are  ill-posed,  in  general,  with  all  three, 
existence,  uniqueness  and  continuous  dependence  on  the  data 
open  to  question.  We  expect  that  these  difficulties  can  only 
be  resolved  when  Y(x)  is  assumed  to  be  non-zero  only  over  a 
finite  domain  interior  to  J5  and  J)  itself  is  a finite  domain. 
Indeed,  because  of  this  latter  constraint,  it  is  expected  that 
the  result  of  Section  II  cannot  be  derived  by  simple  limiting 
arguments  from  (3.12). 
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IV.  LIMITED  APERTURE  OF  OBSERVATIONS 


We  return  now  to  (2.19)  and  assume  that  we  have  only 
a "limited  aperture"  of  observations;  that  is,  we  assume  that 
?(>c)  is  known  only  over  some  finite  domain  as  shovn  in  Figure  4.1. 


I 

I 

1 FIGURE  4.1 

Domain  in  which  ?{jc)  is  known. 

Let  us  define  the  characteristic  function  of  that 

t 

I domain  to  be 


1 , ^9  (jci  is  known 

a(tc\  = • 

0 , ?(<)  is  not  known 

and  denote  by  a(x)  its  Fourier  transform.  Then,  in  fact,  what 
we  really  know  is  the  function  £ (k.) 9 (<) . By  the  Fourier 
convolution  theorem,  this  product  is  the  Fourier  transform  of 
the  convolution  of  y (x)  with  a (x) : 


*(2i)3  J £(<)  ?(<)  exp [ -3k  *x  ] d <3  = J a (x-x ' ) y(x')  dx' 


(4.2) 


■% 
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Here,  the  left  side  is  known  and  this  is  a Fredholm  integral 
equation  of  the  First  Kind  with  the  kernel  a(x)  being  a 
generalization  of  the  cin  x/x  kernel  whicit  would  arise  in  one 
dimension.  Again,  we  note  t-.at,  in  general,  these  integral 
equations  are  known  to  be  ill-posed  with  difficulties  as  to 
all  three  aspects  of  existence,  uniqueness  and  continuous 
dependence  on  the  data.  Nonetheless,  for  this  particular 
equation,  the  second  author  Bojarski  (1973,  3 974)  has  been 
successful  at  obtaining  solutions. 


# 


We  remark  further  that  once  the  idea  of  multiplication 

of  y(k)  is  introduced,  it  may  be  exploited  in  another  way.  In 

particular,  we  note  that  is  is  not  y (x)  that  we  seek  really,  but 

its  derivatives.  Indeed,  it  is  these  functions  that  describe 

the  surface  of  the  scatterer  rather  than  the  scatterer  itself. 

In  Bojarski  (1974),  this  idea  is  exploited  to  develop  a check 

on  the  method.  The  function  p(tc)  is  calculated  from  the  exact 

Mie  series  solution  for  scattering  by  a sphere.  The  inverse 

transform  of  itcj?{£)  is  calculated,  this  function  yielding 

g 

y(x).  Figure  4.2  is  a graph  of  this  function  in  the  (xi,  X2) 
plane. 


The  formulation  of  Section  III  can  also  be  modified 
to  take  account  of  limited  aperture.  To  this  end,  let  us 
suppose  that  3J)  in  Figure  3.1  is  decomposed  into  two  parts 
3pi,  3J>2  as  in  Figure  4.3. 


FIGURE  4.  3 
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On  api  the  data  u and  ^ are  observed,  on  3f>2  they  are  not. 

We  return  to  the  definition  of  0 in  (3.9)  and  decompose  the 
domain  of  integration  into  3Pi  + 3P2;  that  is,  rewrite  (3.2)  as 


0{x,k) 


0)  Q.) 

0 (x,k)  + 0 (x,k) 


(4.3) 


Here 


0<3,<x,k)  - fcs(x\k')  (x-x1)  -g*(*-x')|HS  (x' ,k)|dS' , j = 1,  ?„  (4.4) 


9Dj 


In  0^,  where  we  do  not  know  the  data,  we  replace  ug  and 
by 


3u< 


Ug(x',k)  = & g(x",k)  jjSL  (x"-x' ) - g(x’’-x')  ^j5  (x",k)  j.  dS" 


(4.5) 


and 


***  (x"-x')-^  (x”>zx')  ^ (x\k)|  dS" 


(4.6) 


As  in  Section  III,  we  can  make  the  physical  optics 
approximations  here  and  repeat  for  0*(x,-k).  The  result  is 


0(2)  (x-k)  + fP  (x,-k)  = 0^  (x,k) 


= - IdS'  / dS"fe-  (x-x' 

J J L3n 


) J*  )u  (x" ,k)  g (x"^x‘) 

3^  ) I 


(4.7 


3^2  L&D 


-g*(x-x')  ^ juj  (xM  ,k)  (x"-x')J 


■jS 


a 


Again,  we  can  replace  the  surface  integral  by  a volume  integral 
The  result  is 


(X/k)  = JyW')  K(x",x)  dV" 


(4.8 


Here 


K(x",x)  = / ig (x"-x' ) |2£  (x-xf)  - g* (x-x* ) 3^  (x"-x')|  dS' 
d°2 

Now  the  integral  equation  (3.11)  becomes 

/[*  I?!  (x-x' ) exp  (ik*x‘ ) | + K(x'  /X)l  Y(x')  = 0(^(x,k), 


(4.9 


(4.1 


Equation  (4.10)  is  again  a Fredholm  integral  equation 

of  the  First  Kind  for  6 (x) . It  is  expected  that  the  equation 

becomes  more  and  more  ill-posed  as  8P2  is  the  entire  surface 
u) 

3Pz  = 0 and  the  equation  is  homogeneous.  In  the  opposite 
limit  3Pi  = 3p  and  this  equation  reduces  to  (3.11). 


with  g the  free  space  Green's  function  (2.6).  We  are  justi- 
fied in  restricting  the  domain  of  integration  to  P here 
because  p(x)  is  non-zero  only  inside  of  f>.  Let  us  now 
consider  a volume  integral  similar  to  the  one  used  to  derive 
(5.2)/  but  with  g replaced  by  g*/  namely  we  set 


0(s>  - /[< 


g*  (x-x1 ) Lu(x')  - u(x')  Lg*(x- 


x-x' ) d V' 


(5.3) 


By  using  the  property  that  Lg* (x)  = 3 (x)  , we  find  that 


0(x)  = u(x)  - J g* (x-x' ) p(x' ) dV' , x £ D 


(5.-4) 


By  using  Green's  theorem,  we  also  find  that 


0(x)  = -/ju(x')  |3^  (x-x' ) - g*(x-x')  ^ (x')j 


(5.5) 


Here,  the  integrand  is  known  for  u and  ^ observed  on  3p. 

We  remark  that  a similar  integral  appears  in  the  derivation  of 
(5.2).  However,  by  invoking  the  outgoing  property  for  both  u 
and  g,  one  can  show  that  the  surface  integral  is  zero  in  that 
case;  see,  for  example,  Morse  & Peshbach  (1953)  or  Baker  & 
Copson  (1950).  Thus,  the  process  that  yields  so  simple  a 
result,  (5.2),  for  the  direct  problem,  annihilates  the  data 
that  is  known  for  the  inverse  problem,  while  the  same  process 
applied  to  (5.3)  does  not  do  this. 

We  use  (5.2)  and  (5.5)  in  (5.4)  now  to  obtain 


J 9i 


;-x' ) p (x' ) dV'  = 0 (x)  , 


(5.6) 


As  a check  on  (5.6) , Bojarski  (1974)  has  carried  out 
two  examples.  In  each,  the  function  0(x)  as  defined  by  (5.5) 
is  generated  from  a known  source.  The  check  then  consists  of 
solving  (5.6)  for  p(x)  and  comparing  with  the  known  source. 

The  technique  of  solution  is  a relaxation  iteration  method. 

The  sources  treated  are  a point  source  and  then  two  point 
sources,  one-half  wave  length  apart.  Results  after  8,  16 
and  32  iterations  are  shown  in  Figures  5.2  and  5.3,  respectively. 
The  second  example  strongly  suggests  that  the  method  is  not 
Raleigh  diffraction  limited. 
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RECONSTRUCTED  POINT-SOURCE-PAIR 


VI 
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LIMITED  APERTURE  FOR  THE  INVERSE  SOURCE  PROBLEM 


We  now  suppose  that  u and  are  observed  only  on  a 
portion  of  3p,  namely  on  3Pi  in  Figure  4.3.  As  in  Section  VI , 


we  define  0^  and  0^  by  (4.3)  and  (4.4)  with  0(x,k)  and  u_ 

“ " b 


iu(x' fk)  replaced  by  0(x)  and  u(x' ),  respectively.  In  the 
integral  over  3f>2  we  use  (5.2)  to  obtain  an  expression 
completely  analogous  to  (4.8),  namely  that 


0 (2>  (x) 


x)  = Jp  (x" 


) K(x",x)  dV", 


and  K(x",x)  is  given  by  (4.9).  We  then  obtain  instead  of 
(5.3),  the  modified  integral  equation 


yjgj(x-x')  + K(x',x)j  p(x')  dV'  = 0^(x); 

0^(x)  = j ju(x')||£  (x-x1 ) - g*(x-x')  1 J dS * . 


3p, 
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(6.1) 


(6.2) 


ANALYSIS  OF  INHOMOGENEITY  OF  A MEDIUM 


Let  us  suppose  that  in  an  otherwise  homogeneous 
medium  there  is  a finite  region  of  inhomogeneity. 


FIGURE  7.  1 

Our  objective  is  to  analyze  the  inhomogeneity  from  remotely 
sensed  data.  In  this  section  we  shall  show  how  this  problem 
can  be  reduced  to  the  inverse  source  problem  (5.1)  and  thus 
that  the  inverse  problem  here  reduces  to  the  study  of  equation 
(5.6). 

We  suppose  that  the  wave  (2.1)  (exp  }ik*x|) 
incident  on  the  region  of  inhomogeneity  and  that  the  total 
field  Uj.  + satisfies  the  equation 


(V2  + k2)  (Uj  +us)  = -v(x)  (Uj  + us)  = -p(x) 


Here  v(x)  is  assumed  to  have  finite  support. 


'j€  i,1*  tt$>*  — **4>?  *J*c-s»r  %*  >>v>«®^y  ^ ’ 


In  fact,  Uj  satisfies  the  homogeneous  Helmholtz 
equation  in  the  uniform  medium.  Thus,  (7.1)  reduces  to 


(V2  + k2)  us  = -p(x) 


(7.2) 


3u, 


We  observe  ug  and  j^S  on  3p  and  obtain  the  integral  equation 
(5.6)  for  p.  W!  m p is  known,  then  Ug  can  be  calculated  by 
using  (5.2).  Since  Uj  is  known,  we  can  now  use  the  right- 
hand  equation  in  (7.1)  to  find  v(x). 


In  problems  in  which  one  medium  is  imbedded  in 
another  (tumors,  land  mines,  sea  mines,  buried  objects  in 
general),  v(x)  is  discontinuous  and  (7.1)  is  valid  every- 
where except  on  the  interface  between  the  media.  On  that 
interface,  certain  jump  conditions  must  be  satisfied.  These 
conditions  are  simply  the  integral  form  of  (7.1)  across  the 
interface.  As  a consequence  of  this,  one  can  verify  that 
even  in  this  case,  p(x)  must  satisfy  (5.6). 
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VIII.  THE  TIME-DEPENDENT  INVERSE  SOURCE  PROBLEM 


We  suppose  now  that  u (x,t)  is  a solution  of  the  wave 


equation 


= V2u  - -i-  u = -p(x,t). 


with  U/  U.  and  P assumed  to  be  zero  initially  (sufficiently 

su 

far  in  the  past) . We  assume  that  U and  are  observed  on  3p 
of  Figure  5.1  and  that  P is  confined  to  the  interior  of  p. 

We  shall  derive  an  integral  equation  for  P. 


We  begin  by  observing,  analogous  to  (5.2)  that 


U(x,t)  = 


J f P(x'  »t* 


) G(x-x',t-t')  dV’  dt‘,  t4T. 


Here  G is  the  "causal"  Green's  function.  We  also  introduce 
the  "effectal"  Green's  function  G*: 


G(x, t)  = ?J_t-x/c)  , G*(x,t)  = 6<t  + x/c) 
4irx  ~ 4ttx 


As  in  Section  V,  we  consider  the  function 


IjfeAM-  -*£&&.  *L 


y*  y*  Gj  (x-x'  /t-t*  >p  (x' , t* 


) dV'dt' 


0 J> 


(8.9) 


= 0(x,t,T)  + ${x,t,T),  X Ep,  0<t<T. 


Here 


Gj(x,t)  = G(x,t)  - G*(x,t) 


(8.10) 


If  the  source  is  "shut  off"  after  a finite  time  or 
decays  "sufficiently  rapidly",  then  rp^  $>  = o. 


In  this  case,  we  may  let  T +m  in  (8.9)  and  obtain  the  inverse 
transform  (k  = w/c)  of  (5.6)  as  the  equation  for  the  source  in 
space  time.  When  P does  not  vanish  for  T -*•“  we  are  faced  with 
the  additional  difficulty  of  having  only  a limited  aperture 
(0,  T)  of  observations  in  time.  We  deal  with  this  difficulty 
as  in  previous  sections  by  using  (8.2)  in  (8.8).  The  result 
is  that  P(x,  t)  must  satisfy  the  integral  equation 


T 

fj 


\x (x-x'  ,t-t' ) + K(x,x',t,t';T)j  P(x',t')  dV'dt’ 
= 0(x,t,T),  x e«p,  o<t<T 


(8.11) 


and 


K(x,x',t,t';T)  = y|c(x'-x",T-t')  _jj_  G*(x-x",t-T) 

P (8.12) 

-G*  (x-x" , t-T)  ^1-  G (x'  -x"  ,T-t ' ) | dV" . 


One  can  certainly  modify  this  equation  further  to  take 
account  of  a limited  aperture  of  observations  in  space.  However, 
this  modification  is  carried  out  exactly  as  in  Section  IV  and 
VI  and  we  shall  not  do  so  here. 
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